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ABSTRACT 


Scour About A Single, Cylindrical 
Pile Due to Random Waves and Current. (December 1984) 
Robert W. Eadie IV, B.S., United States Naval Academy 


Chairman of Advisory Committee: Dr. J. B. Herbich 


An empirical approach was utilized in this investiga- 
tion to study the scouring effects that various currents 
and random waves have on a single cylindrical pile. Re- 
sults from laboratory experiments are presented in graphi- 
cal form depicting relationships between pertinent dimen- 
sionless parameters which were developed through a dimen- 
sional analysis. These results were then compared to re- 
sults compiled in previous studies that investigated scour 
due to currents and monochromatic waves in order to deter- 
mine if the irregularity of the wave train is significant 


in the study of scour. 
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PART I 


INTRODUCTION 


Eve COWsErUcCtIon Of Offshore Structures is an ex- 
tremely expensive process. A portion of these construc- 
tion costs are due to large safety factors used during 
the design process. These safety factors which may lead 
to overdesign and, therefore, unnecessary expense are the 
result of the engineer's inability to accurately predict 
some of the forces on the structure, as well as some of 
the natural phenomena that will occur once the structure 
is in place. One of these natural phenomena which is be- 
coming an ever increasing problem is scour caused by waves 
and currents around offshore structures. Design factors 
for offshore structures depend on the location of the ocean 
bottom. The scour of a non-cohesive ocean floor can lower 
the bottom line significantly, and by doing so can dras- 
tically change these design factors. 

The purpose of this investigation was to study what 
effect various currents and random waves have on a singu- 
lar, cylindrical pile. These results were then compared 
to results compiled in previous studies that investigated 
scour due to currents and monochromatic waves to determine 
if the irregularity of the waves is significant. An em- 


pirical approach was used in this study. 





Experiments were conducted in a two-dimensional wave 
tank. The pile used in this study had a diameter of 1.5 
inches. Two water depths, two currents, one sediment size 
and four random wave spectra were utilized. Using data 
obtained from the experiments, an attempt was made to de- 
scribe scour in terms of relevant dimensionless and dimen- 
sional parameters. By doing so, it is intended that the 
reader will obtain a better understanding of the scouring 


process. 





PART II 
REVIEW OF LITERATURE 


Incipient Motion 


The scour phenomenon begins with the incipient motion 
of sediment particles. Eagleson and Dean (7) define in- 
cipient sediment motion to be "an instantaneous condition 
reached when the resultant of all the active forces on 
the particle intersects the line connecting the bed parti- 
cle contact points." There are three types of active 
forces acting on a sediment particle. Two of these forces, 
the drag and lift forces, are due to the fluid motion. 
These forces are in turn balanced by the sediment parti- 
cle mass. These forces are depicted in Figure 1 (1). 

A thorough understanding of these forces is required 
in order to comprehend fully the processes of incipient 
motion, bed-load transport and bed form development; all 
which depend on the collective response of initially sta- 
tionary bed particles to fluid forces. The following dis- 
cussion of incipient motion is for a unidirectional flow. 

The variation of drag forces with flow condition has 
been well documented, as has the lift force in hydrodyna- 
mically-rough flow conditions. However, investigations 
of lift forces in hydrodynamically smooth and transitional 


flow have yielded conflicting results. 
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Prior to 1967, most investigations measuring fluid 
lift forces were conducted in the hydrodynamic rough flow 


region characterized by: 





in which R, represents the grain Reynolds number; V,, 
shear velocity; D = grain diameter; and v = kinematic vis- 
cosity of the fluid. All measurements under these hydro- 
dynamic conditions indicated Sosiieae dit fae. 

Coleman (4) was the first to study lift forces in 
the hydrodynamically smooth region (R, < 5) and transition- 
al flows (5 < R, < 70). His work reported measured drag 
forces and inferred lift forces from observing the inci- 
pient motion of spheres having various densities. From 
the results, it was suggested that for a sphere resting on 
top of a closely packed bed of spheres of equal diameter, 
the lift force decreased through zero ann became negative 
(bedward-directed) as R, decreased through a value of ap- 
proximately lL5. 

Davies (6) prepared an experiment whose sole purpose 
was to define the flow conditions under which negative 
lift forces occur. He found that significant negative 
lift can be expected to occur only if the following con- 
ditions are met: 1) The value of R, < 5 and 2) The test 
sphere is situated such that a significant volume rate of 


flow can occur beneath the sphere. Ina situation in 





which considerable Fluctuation of bed shear stress occurs, 
the onset of significant negative lift with a decrease in 
R,., will occur gradually and unsteadily as the effect of 
turbulence in the area of the bed particle becomes less 
significant, until hydrodynamically smooth flow is es- 
tablished at R, < 5 and the lift force becomes predominant- 
ly negative. 

The existence of varying types of lift forces is im- 
portant with regard to bedform development. If the par- 
Erelesis removed from its position of rest at Ry < 5, the 
negative lift forces acting on the particle will tend to 
prevent the particle from leaving the immediate vicinity 
of the bed, and, therefore, it will move along the bed as 
a bed load. On the other hand, if the particle is removed 
at R, >> 5, it will experience a positive lift force that 
will tend to impart an upward acceleration and thereby 
lift the particle into the main turbulent flow, where it 
may or may not become suspended. 

There is considerable evidence to suggest that a rip- 
pled bedform exists only at R, < 5, and that a dune bedform 
exists on its own only at R, > 24. It is therefore sug- 
gested that: 1) The flow conditions that are necessary 
for existence of the negative lift forces on a sediment 
particle correspond closely to those required for a plane 
or rippled bedform. 2) Positive lift force conditions 


correspond to a dune bedform. 





The condition of incipient motion for a cohesionless 
soil particle is described in terms of the forces acting 
on the particle by: 
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where F and FA are the forces parallel and normal to the 


angle of repose 0. In our case, F, and FA are the resultant 


IE 
forces of the three active forces mentioned earlier, the 
hydrodynamic drag force Fa the lift force FL and the 
submerged weight W (gravity force). In accordance with 
equation (1), the condition of incipient motion under the 
influence of these three forces becomes: 


Wsina + F 
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where the angle a is the inclination of the bed from the 
horizontal at which the incipient motion occurs. The above 


situation is illustrated in Figure 2. 
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Fig. 2. ~ Force diagram on particles in a cohesionless 
loose bed (9) 
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The submerged unit weight of the particle is expressed 


as: 
W= K.,(p. - p)gd- 
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By introducing equations (3), (4) and (5) into equation 


(2), the result becomes, 


(U,*) Ong 2K, (tan © cos a- sina) 
— aa aCe ean oe (6) 
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where (U,)er is the critical bottom velocity at which ac- 


emt elo equation (1) incipient motion occurs. 
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where 


T eritical shear stress 


O 

0 = fluid density 
The value of the right hand side of equation (6) de- 
pends on A) The particles: their size, shape, size dis- 


tribution, uniformity, .. .; B) The flow dynamics, which 


determine the values of Ch and C, C) the channel slope 





which in many cases amounts to cos a * 1; and D) The angle 
of repose which depends on the particle properties (9). 

Once Ehercritical@eliow velocity#is reachedesthe cri- 
tical shear stress will also be obtained as per equation 
(7). Most data on critical shear stress for non-cohesive 
sediments have been observed in flume experiments. These 
experiments have shown that the motion of the sediment 
particles at the bed for uni-directional flow is highly 
unsteady and nonuniformly distributed over the entire bed 
area. Near critical conditions, the particle motion in 
any small area of the bed occurs in gusts whose incidence 
increases as the mean shear stress increases. When ob- 
serving a large section of a sediment bed as the shear 
stress approaches the critical value, the occurrance of 
these gusts of particle motion seem to be random in both 
time and space. 

Two of the best known and widely accepted investiga- 
tions into the concept of critical shear stress are those 


of Shields and White (17). Shields developed the relation: 
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) is shown in Figure 3. 

As discussed earlier under Davies experiment, under 
flow conditions in which R, < 5, the particle will have 
negative lift and will therefore be less likely to move. 
This corresponds to the Shields curve that indicates that 
when R, decreases below about 10, particles increasingly 
resist the dimensionless shearing stress of a flow. 

White (1/7) determined the critical shear stress for 
several sediments under conditions of developed laminar 
flow, laminar boundary layer flow and turbulent boundary 
layer flow. White concluded that the true critical shear 
stress necessary to move a sediment particle has a fixed 
value and that this value is given by experiments dealing 
with laminar flow. White's (24) equation for we for sedi- 


memes in a horizontal bed is 


tT, = 0.18 Cope E20 Osa ee ee (9) 


where 


critical shear stress (dynes/cm*) 


j 
It 


the angle of repose of the sediment submerged 


in the fluid 
and in which the constant (0.18) is obtained from experi- 


ments with laminar flow. 
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Fig. 3. - Shields diagram with White data added (24) 
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Yang (28) expressed some reservations on using Shields' 
diagram as the criterion for incipient motion. Using basic 
concepts in fluid mechanics, the boundary layer theory 
and laboratory and field data from various investigators, 
Yang plotted the relationship between dimensionless criti- 
cal velocity and shear velocity Reynolds number at inci- 
pient motion (Figure 4). From this curve, the following 


equations of incipient motion were obtained: 











V Diesel 
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log (=) - 0.06 
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where 

V = average critical velocity at incipient motion 


CR 
W = terminal fall velocity, which is reached when 


the particle drag force and the submerged unit 
weight of the particle are balanced. 
d = particle diameter 
“= kinematic Viscosity 
U. = shear velocity (U,. = ae oe 5 Lic Cenisat yy 
t = average shear stress) 


U,d 


= shear velocity Reynolds number 





Many experiments have shown that for any given bed 
material, there is no single value of bed shear stress or 


velocity, below which not a single particle will move and 
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above which all the particles having the same size will 
move. The forces acting on the particle are not constant, 
but instead fluctuate about some average value, a phenomenon 
closely associated with turbulence. Also, the forces which 
move the particle are dependent on the particle position. 
Therefore, as mentioned earlier, the movement of bed ma- 
terial particles is very unsteady, and non-uniformly dis- 
tributed over the entire bed (17). Einstein (8) has sug- 
gested that the process of bed material movement is sta- 
tistical in nature. Because of this, it appears that no 
truly critical condition for particle movement exists. 

This is further exemplified in Figures 5 and 6 which show 
the variance in values for the critical shear stress and 


critical water velocity for quartz sediment. 
Scour 


Fluid Mechanics of Scour 
Local scour is defined as "the abrupt decrease in 
bed elevation near a pier due to erosion of bed material 
by the local flow structure induced by the pier." (22) 
One type of classification of the type of scour is the 
amount of sediment transport into and out of the scour hole. 
Laursen (12) suggested that the scour process could 


be represented by the relationship: 
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in which dQ. /dt pele babe sor loeal scour) inevo lume sper 
unit time; ae the capacity of the flow to transport 
sediment out of a scour hole in volume per unit time; and 
dg2 = the rate at which the undisturbed flow supplies sedi- 
ment to the scour hole. 

From this expression, three cases of scour may occur. 
They include: 1) no scour, in which case Pee Sea 5. = G: 

2) clear water scour, where Sl se O and Q ~ G59 << 4613 and 

3) scour with continuous sediment motion, where eae Ge ON 
Abrupt changes in bed elevation near a pile occur only in 
cases (2) and (3). These two cases of scour are depicted 

in Figure /7. 

The main flow feature responsible for scour whenever a 
slowly-varying boundary layer current is present, results 
from secondary flow. The passage of a steady flow boundary 
layer along the surface past the cylinder causes a stagna- 
tion pressure gradient. Because the stagnation gradient 
decreases ‘toward the sediment bed, a secondary flow is 
generated in the direction of the bed. This secondary flow 
has two effects on the shear or boundary layer flow. The 
principle effect is the formation of a horseshoe vortex 
that wraps itself about the cylinder and trails off down- 
stream. The horseshoe vortex shown in Figure 8 is depicted 
as a spiral curve extending around the cylinder base (15). 

Considerable laboratory data are available from scale 


model studies; however, most of these studies investigate 
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the horseshoe vortex in a river environment. Even though 
these results provide some insight into the characteris- 
tics of flow and the basic mechanism of generation of the 
horseshoe vortex, it is impossible to apply these river 
model results directly to an ocean environment. A major 
deficiency being the lack of appropriate dynamic and kine- 
matic scaling laws. Shen et al. (21) have shown that the 
size and the strength of the horseshoe vortex can be re- 
lated to a Reynolds number using the characteristic length 
of the cylinder in the flow field. This relationship is 
as follows: 


Zau 


Vortex strength = “= f(-——") .... 2... 2. s (34) 





where 

i esorattonal velocity of vOrtEex core 

A = area of vortex core 

a = half width of pile 

Be = undisturbed flow used in the potential flow sense 

vy = kinematic viscosity of fluid 

Also, Shen et al. (21) sue rested that for cylinders 
with a height greater than the thickness of the bottom 
boundary layer, the area affected by the horseshoe vortex 
extends approximately one cylinder diameter upstream om 
@me Obstruction. 


Hjorth (10) indicated that the fluid shear stress 


striking the bottom below the horseshoe vortex system are 





increased by a factor of eleven when compared with the 
shear stress values associated with the undisturbed free- 
stream. This increase in bottom shear stress causes the 
bottom sediments to scour. With scour occurring, the bot- 
tom elevation is altered, which in turn brings about a 
change in the flow in the horseshoe vortex. 

Another effect of the secondary flow about a cylindri- 
cal obstruction is the downward deflection of the external 
Flow streamlines about the cylinder. This downward deflec- 
tion of streamlines is shown on Figure 8, the secondary 
flow is depicted at the stagnation line on this same fig- 
ure. The convective effect of the streamlines has little 
effect on the horseshoe vortex, except for a slight strength- 
ening of the rotational vortex. 

Pecitculareeylinder blocks the approach Elow, result- 
ing in convective acceleration. The acceleration causes 
an increase in the bottom shear stress by a factor of up 
to four over the bottom shear stress of an unobstructed 
Flow. 

Flow patterns in the cylinder wake also influence 
scour. There is a net downward flow in the wake of a 
cylinder resting on a plane surface. This is due to the 
fact that the wake region is a low pressure region. The 
wake in the high velocity region is at a higher pressure 
than the wake in the low velocity region near the sedi- 


ment bed. This difference in wake pressure causes a flow 





to form in the wake acting downwards toward the bed. If 
the effects of the horseshoe vortex are included, this 
downward flow is enhanced (15). 

In an attempt to get a better understanding of the 
mechanisms of scour, the subject will be dealt with in 
three parts. First, the mechanics of scour in a steady 
current alone will be discussed, followed by the mechanics 
of scour due to oscillatory flow, and, lastly, the mechanics 
of scour caused by the interaction of waves and currents. 
Mechanics of Scour in Steady Currents 

Repeated observations indicate that the scale and 
shape of the scour hole closely resemble the scale and 
shape of the horseshoe vortex. The scour is caused by 
the large bottom fluid shear stress gradients occurring 
beneath the horseshoe vortex. The initial vortex, before 
the onset of scour, exerts bottom fluid shear stresses 
that are approximately one order of magnitude greater than 
the shear stresses in the surrounding undisturbed flow. 
Therefore, it is possible for scour to occur even though 
sediment is not eroded and transported into surrounding 
regions. If general bedload transport is occurring, then 
the elevated bottom shear stresses beneath the horseshoe 
vortex produce more rapid sediment transport and bed ero- 
sion. 

As the erosion of sediment continues, the intensity 


of the horseshoe vortex decreases as it spreads out into 





the scour cavity. As the vortex intensity decreases, the 
fluid shear stress acting on the sediments in the bottom 
of the scour cavity are also decreased. Once equilibrium 
is reached, the fluid shear stress acting upstream at the 
bottom of the scour cavity is equal to the fluid shear 
stress acting downstream on the undisturbed portion of 
the bed. The fluid shear stress acting on the upstream 
face of the scour hole is equal to the fluid shear stress 
acting on the undisturbed area of the sediment bed, plus 
an increment approximately equal to the fluid shear stress 
required to entrain the sediment. Once this distribution 
of shear stresses has been achieved, the scour hole has 
reached its equilibrium size and shape (Figure 9) and the 
rate of transport of sediment from the scour hole equals 
the rate of supply (15). 

Mechanics of Scour for Oscillating Flow 

Observations of scour development from oscillating 
flow are given by Rance (19) and Palmer (18). The exact 
sequence of scour hole development depends on the diameter 
of the cylinder, the intensity of the oscillating flow, 
and the sediment transport on the surrounding bed. 

Rance states that for a slender body (i.e., the dia- 
meter of the body is small when compared with the wave- 
length) under the action of waves, the shape of the scour 
hole is dependent on the amplitude of the wave motion at 


the bed. When the amplitude is very small when compared 
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Fig. 9. - Steady current scour cavity clos 





to pile diameter, two discrete scour holes are formed on 
the forward part of the cylinder at 45° to each side of 

the centerline. For this case, eddy shedding has no ef- 
fect on sediment transport. When the bed semi-orbit length 
is larger and becomes approximately one fifth of the pile 
diameter, there is a tendency for flow separation to occur. 
Accretion and erosion are localized and four scour holes 
versus two are formed. These scour holes have a maximum 
depth of about 10% of the pile diameter. At larger ratios 
of semi-orbit length divided by pile diameter, of about 
unity and above, eddy formation is pronounced and a scour 
hole develops completely and uniformly about the pile with 
a maximum depth of 40% of the pile diameter. 

Palmer noted that the scour hole initiates at the 
sides of the cylinder and, in many cases, some sediment 
deposition was found to occur against the upstream and down- 
stream sides of the cylinder. As the scour hole continues 
to develop, these deposits are removed and a majority of 
the sediment deposition shifts to a rim outside the scour 
hole. Once fully developed, the scour hole was found to 
be a radially symmetric frustrum of an inverted cone hav- 
ing side slope at the angle of repose of the sediment. 
Palmer's data indicated that the depth oF the equilibrium 
scour hole to the pile diameter ranged from between ().5 


and 0.75, a value higher than that previously stated by 


Rance. 





Horeosca lating Llow past a cylindrical obstruction, 
the horseshoe vortex is insignificant, and, therefore, 
scour is initiated from the bottom fluid shear stress 
caused by the acceleration of primary flow past the ob- 
struction. This results from a bottom fluid shear stress 
gradient that is positive from upstream of the cylinder 
to the point of maximum velocity, and negative from that 
point to the downstream limit of undisturbed flow (Figure 
10). If the maximum bottom fluid shear stresses are great- 
er than the sediment entrainment threshold, erosion will 
occur from upstream of the cylinder up to the point of 
maximum velocity, and deposition will occur downstream of 
ears point. 

In order for an equilibrium scour hole to develop 
around the entire circumference of the cylinder, an addi- 
tional flow condition must develop. The initial scour hole 
produces bottom slopes close to the pile which must even- 
tually cause flow separation in the thin orbital boundary 
layer. Reverse eddies in the separated flow erode the 
initial side slope of the scour hole until it obtains a 
value of the angle of repose of the sediment. With the 
steepening of the sideslopes, the zone of separated bottom 
boundary layer flow spreads out around the pile eventually 
causing the scour hole to obtain the equilibrium shape 
shown in Figure 11. Therefore, the equilibrium scour hole 


for oscillatory flow is related to boundary layer separation 
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Fig. 10. - Sediment transport due to oscillatory flow (15) 
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and not dependent on the nature of secondary flow as was 
the case for scour due to steady currents (15). 
Mechanics of Scour Due to Waves and Currents 

In this study, the scour that developed through the 
combined action of waves and currents was investigated. 
This is the type of scour that occurs in almost all ocean 
environments. Ninomiya et al. (16) performed scale model 
experiments using the combined action of waves and currents. 
They found the shape of the scour hole to be similar in 
shape to that developed by a steady current without waves, 
and that the rate of scour hole development was consider- 
ably faster in the case of combined wave and current action 
when compared to steady current. 

As mentioned earlier, if the steady flow velocity is 
stronger than the oscillating flow velocities, a horseshoe 
vortex will form about the base of the pile. The intensity 
of this vortex varies according to differences in the 
square of the-combined steady and oscillating flow velo- 
cities. The size and shape of the equilibrium scour hole 
should not be as large as that caused by a steady current 
whose mean velocity is equal to the combined velocity of 
the wave orbital and mean current condition. This is due 
to the fact that velocity profile of the combined veloci- 
ties will be flatter than that achieved in a simple 
steady boundary layer flow. On the other hand, the equi- 
librium scour hole depth should be greater than that caused 


by a steady boundary layer having no wave El@lealopoe 





If a strong current and a weak wave orbital flow 
exists, the development of the scour hole will depend 
primarily on a horseshoe vortex with a pulsing intensity. 
The scour hole development will be faster when compared 
to a steady current alone due to the additional bottom 
fluid shear stress caused by the wave orbital boundary 
layer and because the vortex formed from the steady por- 
tion of the flow can more easily transport the wave- 
agitated sediment. 

In the case of a strong orbital flow and a weak cur- 
rent, the scour hole shape will be similar to that caused 
by a steady current. Once again, the scour will develop 
more rapidly than it would under the influence of a steady 
current alone due to the additional bottom fluid stresses 
caused by the waves orbital motion. In this case, the 
equilibrium scour depth can be less than that associated 
with a steady current that is equal to the sum of the weak 
current and orbital velocity. 

The situation becomes difficult to evaluate when the 
magnitude of the oscillatory motion and the steady current 
are nearly equal. The scour hole size and shape will nor- 
mally depend on the horseshoe vortex associated with a 
steady current. It is not possible to generalize on the 
scale and intensity of the pulsing vortex when compared 


to the steady current situation. 





The size and depth of scour produced by the combination 
of waves and currents acting on a pile that is large in 
comparison to the oscillation amplitude will be drastical- 
ly smaller than that predicted for a small pile, except 
when the steady current is greater than the maximum wave 
orbital velocity. This is caused by the fact that the 
horseshoe vortex does not have time to wrap completely 
around the pile before the flow reverses due to the wave 
Spoiled! velocity (15). 

Previous Studies on Scour 

Very little experimental work has been done on scour 
caused by the combined effect of oscillatory wave motion 
and current; however, there have been numerous studies 
dealing with scour in open channel flow. Because the 
scouring forces are similar for oscillatory flow as for 
steady state open channel flow, that is, hydrodynamic in 
nature, the results found in experiments in open channel 
flow can be applied with reservation to the combined oscil- 
latory and steady current flow. 

Some experimental results are presented by Shen et 
al. (22) as shown in figures 12-14. These experiments 
employed a steady current. Figure 12 indicates that the 
depth of scour is a function of the pier Reynolds number 


Gina. 
R = UD/v 


where 
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G 
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horizontal free stream velocity 


D 


pile diameter 

Vea tnematiec viscosity 
It was also shown that the scour depth is also dependent 
on the pile diameter (Figure 13) and sediment size (Figure 
14), and this explains why the Reynolds number alone is 
insufficient to collapse all of the data shown in Figure 
12 onto a single line. 

Carstens (3) performed studies of scour associated 
with various types of obstacles. He concluded that the 
rate of scour induced by the obstacle is a function of 
the sediment number (No), sediment geometry, as well as 
the ratio of the scour depth to obstacle size. It should 
be noted that the ''sediment number" is a densimetric Froude 


number. 


N. = U/ (S.-1l)gd 


S 
where 
U = free stream velocity 
S. = sediment specific gravity 
g = acceleration of gravity 


d = mean diameter of sediment particle 

Coleman (5) used the principle that in any problem 
of flow around an immersed body, where the dependent vari- 
able is a pressure drop, the most general formulation of 


the problem is in terms of an Euler number and either a 








Reynolds number, a Froude number or both, with the dependent 
variable appearing in the Euler number. He derived a rela- 
tion between the scour Euler number (u/ 2g4_) and the pier 
Froude number (U/ VgD). 


be U .9/10 
aay QO. 6 (—- s s s e e e e s s e e e e eo e e e 6 
Tea = (35) 


From the above relation and using a previously derived 
depth of scour relation using the pier Froude number, an 
expression was obtained to predict the depth of scour for 


steady flow. 


9/10 Ly AL, 


d_ = 1.49 D (U*/2g) (36) 
where 

d. = depth of scour (ft) 

D = pile diameter 

U = free stream velocity 


g = acceleration of gravity 

Wells and Sorensen (27) investigated scour around a 
single circular pile due to oscillatory wave motion. Fig- 
ure 15 shows that a relationship exists between the rela- 
tive scour and the relative water depth for various sedi- 
ments. In addition, their results indicated that the ra- 
tio of ultimate scour to wave height appeared to be di- 
rectly related to the sediment number , No. (Figure 16) and 
to the pile Reynolds number, Np», (Figure 17). In Figures 


Pemenoueh 17, the notation is defined as follows: 
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S. = ultimate significant scour depth 


d = mean grain diameter 


Sana no, tL: d= .62 mm 
sand no. 2: d= ,.325 mm 
band no. 57d = .30.mm 


h = still water depth 

g = acceleration of gravity 

T = wave period 

H = wave height 

U = maximum free stream velocity 
So = Specific gravity of sand 


D = pile diameter 
v = kinematic viscosity of water 
More recent scour studies include the work of 
Ravdkivi and Ettema (20), Jain and Fischer (11), Armbrust 
(2), and Wang (26). The work of Armbrust and Wang dealt 
with the scour associated with the combined effect of waves 
and currents while the remaining studies dealt only with 
a steady current. 
Raudkivi and Ettema (20) performed laboratory experi- 
ments on clear water scour of cohesionless bed sediment 
Sie mandrical piles. They found that; 
1. As the geometric standard deviation of the par- 
ticle size distribution increases then the equi- 


librium scour depth decreases. 





2. The equilibrium depth of scour decreases with 
the decreasing relative size of pier diameter 
to mean bed particle size for values less than 
abomeeZO LOZ). 

3. The maximum equilibrium depth of scour is a 
function of the relative depth h/D and the rela- 
tive pier size D/Dso- 

4. The equilibrium depth of scour decreases at a 
greater rate with decreasing flow depth for 
smaller values of the relative flow depth h/D. 

Jain and Fischer (11) performed scour experiments 

for steady currents at high velocities, i.e., for flow 

regimes well beyond the onset of general sediment trans- 

port. From these studies, they concluded that after de- 
creasing slightly upon the start of sediment transport, 
scour depth does increase again with increasing velocities. 

This conflicts with the consensus of previous experiments 

which concluded that once the conditions of general sedi- 

ment motion have been established in a stream, there will 
be no further increase of scour depth with velocity, be- 
cause the rate at which sediment enters the scour region 

is equal to the rate of sediment removal by the scouring 

process. 

Armbrust (2) found that in general, for a given cur- 

rent and relative depth (h/L), scour increases as the 


wave steepness (H/L) increases (Figure 18). Figure 18 
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also shows that the ultimate scour depth increases as the 
relative depth decreases. Additionally, in general, for 

a fixed relative depth, Figure 19 reveals an increase of 
scour depth with an increase of relative wave height (H/h). 


In Figures 18 and 19, the notation is defined as follows: 


S 


a ultimate scour depth 


D = pile diameter 
Deo = mean sediment diameter 
H = wave height 
L = wave length 
Ui 2 = free stream velocities 
h = still water depth 
Wang (26) defined ultimate scour (SD) to include 
"seneral scour" and net "local scour'’ and, therefore, ul- 
timate scour is equal to the total scour that occurs. Wang 


Stated that: 


5 = £(V__,h,o,H,L,T,do9,9,,D,v,t) 


u 
where 
Ue = combined wave and current velocity 
h = water depth 
09 = water density 
H = wave height with current 
L = wave length with current 
T = wave period with current 
De— pile diameter 
d = mean sediment diameter 


50 
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density of sediment 


12, 
ll 


S 


v kinematic viscosity of water 


t time required to produce scour 


Using dimensional analysis, Wang produced the relation: 


S 
Ay ; : 5 = 
= f(N, H/L UL N. Npp) £(B) 
where 
oe wt BU; No @NSD 
voH*L[V + (1/T - V/L) i" D 
aa” 4 


ae eee 50 


Figure 20 is a plot of the relationship of relative scour 


depth, S/h, and the parameter 8B. 
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THEORETICAL ASPECTS 


Wave-Current Interaction 


Hoften and Karaki (23) performed a laboratory study 
to investigate the interaction of waves and currents as 
mechanically generated, monochromatic waves propagated 
on a turbulent open channel flow. 

They found that waves propagating on an open channel 
flow are shown to be modified by a combination of the wave 
interaction with the shear gradient and the existing tur- 
bulence. The effect was to remove energy from the waves 
causing an increase in wave energy dissipation. 

The flow itself is also changed by the waves. The 
induced orbital wave velocities distort the mean velocity 
profile increasing the gradient close to the sediment bed. 
Turbulence energy is formed near the surface by the in- 
teraction of the wave-induced Reynolds stresses and the 
fluctuating velocity gradients. The energy diffuses down- 
ward where it is dissipated on the sediment bed in the 
form of an apparent increase in boundary shear stress, an 
important factor in the scouring process. 

In order to evaluate the combined velocity of a wave 
propagating in the same direction as a uniform current, 
the Airy wave theory will be used. By doing so, if a 


uniform current is assumed, the velocity potential of the 
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uniform current ($4) and the velocity potential due to 
wave (> 5) can be added directly to obtain the velocity 
potential of the combined wave current condition (¢) be- 
cause both velocity potentials are linear solutions of the 
Laplace equation. The parameters of wave period (T), wave 
height (H) and wave length (L) are shown in Figure (21), 
and are properties of a wave propagating in still water. 
Wave velocities can be calculated using these parameters 
with the water depth for a no current condition. There- 
fore, the velocity relation for combined wave-current in- 
teraction becomes, 


¢ = oy = 5 5 ee Ae Ce ot ee ea | A ee trans ly.) 
where 


oy Ux ee Sere ee re ce CES 
and the derivation of 5 is outlined in the following pages. 
This procedure closely follows Newman (14). 

The fluid velocity under wave action is expressed by 
the gradient of the velocity potential 4, and the effects 
of the free surface must be expressed in terms of appropri- 
ate boundary condition on this surface. 

The exact kinematic boundary condition can be derived 
most readily by requiring that the substantial derivative 


of the quantity z - n vanish on the free surface. The 


result of this condition is that on 2 =n, 


_ D er Mie a Ae eee 1} 
O- FF (z - n) == ned dt x JX dy sy ° . . ° ® ° (14) 
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For the linearized kinematic boundary condition, the second 
order terms can be dropped, which gives the relation: 


en _ 9a 
lp) Oe 2 a a a i aa. 


This approximate boundary condition states that the verti- 
cal velocities at the free surface and the fluid particles 
are equal, disregarding small departures of the surface 
from the horizontal orientation. 

The dynamic boundary condition is obtained from 
Bernoulli's equation. The exact condition to be satis- 


fied on the free surface is: 


1 3 
2 =€ s P.) = =o + eZ ma ee — (ee (1G) 


Substituting the free surface elevation n for Z and solv- 


ing for n the expression becomes, 
= d¢ 
n = -l/g et leon: Vo OIE ea eo, we Se ee ALT) 


By neglecting the last term, which is of the second order, 


the linearized equation for the free surface elevation is 
a chy 
n l/g Se ce an A I eg (18) 


Pimede surtace, 2 — 0, the Vinearized dynamic boun= 
dary condition (18) can be differentiated with respect to 
time and added to the linearized kinematic boundary con- 


dition (15), to form a single boundary condition for the 


velocity potential: 





SJ 1 


2 
dg 9 oy = 
se tS Se 7 OM z= 0... ... ee. ee ee. 9) 


For plane progressive waves with finite-depth effects, 


the free surface elevation will be of the general form, 
fee) = a cos (kx - wt +e) .....,...2..6. (20) 


in which the positive x - axis is in the direction of wave 
propagation, a is the wave amplitude, and « is some arbi- 
trary phase angle that can set equal to zero by the suit- 
able choice of the origin -x = 0. The parameter k is the 


wave number, in which, 
ee 27/7. 


where 

L = wavelength 

The solution to this problem is written in terms of 
a two dimensional velocity potential $(x,z,t), which must 
satisfy. the free-surface condition (19), Laplace's equa- 
tion, the bottom boundary condition 346/3z = 0 at z = -d 
(for finite depth), and must produce the wave elevation 
(20) for equation (18). The potential must be sinusoidal 


due to (20) and the solution should be in the form: 


mee ze ty) = Re(Z(z)e 1 a5 iwt, ee eens 21D) 


From Laplace's equation, Z must satisfy the ordinary dif- 
ferential equation: 


2 

daz 2, | 22 

nee le uz (Cans NU eS 0 er a 722) 
Zz 








The most general solution to Equaiceron (72) inet cams 


of exponential functions is 


emcee - + pa? Os act (OS ) 


where 

C = constant 

D = constant 

First assume the fluid depth is infinite or 
-e < z< 0. To avoid an unbounded motion deep below the 
free surface, the constant D in equation (23) must vanish 


leaving 
ee OD 


If equation (24) is substituted into equation (21), the 


velocity potential becomes, 


6 = Re[CeX@-tketiwt) 88) 


By setting z = 0, using equation (18) to calculate the 
linearized wave elevation, and comparing the result with 
the required plane wave in equation (20), the constant 


C has the value of 
etc, Gs . Cs Cs Ce Ce te «= 26) 


For a fluid having a constant depth d, with the bot- 
tom being an impermeable plane, as stated earlier, the 


bottom boundary condition becomes, 





a> _ = 
— CRC. teers ee OR ek em es s 627) 


By looking at the general solution (equation (23)), it is 
found that both exponential functions must be retained, 
with constants C and D chosen such that equation (27) will 
be satisfied. 

For simplification, the two independent solutions 
e es are replaced by linear combinations of the hyperbolic 


Sine and cosine function. For a fluid having finite depth, 


the solution becomes, 


cosh kk (Z@ +d) 00 0 88) 


Z(z) = © cosh 


The constant C is selected such that equation (26) remains 
valid. When substituting (28) into (21) and using equa- 


tion (26) the velocity potential for a finite depth fluid 


becomes, 
p, = 82 cosh ee OR Gc a .. 


Substituting equation (29) and (13) into equation 


(12), the combined wave-current velocity potential becomes: 


cosh k(z + d) 


aoe SH Keane = © sas + pO) 


= Ba 
o Ux + 


Differentiating equation (30), the result becomes: 


~y 4 gak cosh kK(Z +4) cos xe - wt) . 2. . (32) 
& W cos 


qt Q 
% [o> 


U = 


where 





U = combined wave-current velocity 

It can, therefore, be concluded that the combined 
wave-current velocity (U) can be obtained by simply add- 
ieecme: Value for the uniform current velocity to the 
wave velocity in still water found using the parameters 
stated earlier. 

It is appropriate at this time to mention several 
points. Firstly, this relationship is invalid for non- 
linear boundary conditions. Also, since the wave pro- 
pagates into the flow, there will be modifications to 
both the waves and currents by interaction, therefore, 
the simple addition of current velocity and wave velocity 
is only approximate. Finally, linear wave theory is only 
valid for small amplitude waves, and many cases of wave- 
current interaction exist beyond the accepted range of 
Maiaiey for this theory. 

Wang (26) calculated the combined wave-current velo- 
city using the procedure outlined above and compared the 
calculated values with measured values. He derived a 


combined wave-current velocity relation given by: 


IL 27 
U = Un ate tH(1/T = U/L) sinh 27 d cosh[>— (z ate d) (32) 


L 
where 
U = combined wave-current velocity 
U. = uniform current velocity 


By assuming that the distribution of current is uni- 


form with depth, the depth-averaged velocity was obtained 





from the integral of equation (32). When comparing the 
calculated and measured combined wave-current velocities, 
a correlation coefficient of 0.939 was obtained, thus in- 


dicating this method for calculating currents is reasonable. 


Dimensional Analysis 


There are many variables encountered when studying 
scour. The following list is believed to be the important 
variables influencing scour depth at a cylindrical pile 


gue to cusrent and random wave motion. 


Description Symbol Units Dimensions 

Water depth h ft le 

Significant wave Hi /3 inte L 

height 

Kinematic viscosity v foie L?/T 

of water 

Pile diameter D faite L 

Mean diameter D IIE L 
; 50 

of sediment. 

Density of water p slugs/ft~ M/L> 

Sediment/fluid (0 5 - 2) slugs/ft> M/L? 

density difference 

Acceleration due g Beeece L/T? 

to gravity 

Combined wave- U ft/sec lig JE 

current velocity 

Ultimate scour Sq rte L 


depth 


Since random waves were used in this study, and they 


are defined by the significant wave height (Hy 73)> wave 





a0 


period and wave length have been omitted from the above 
iesteOL variables. 

The dimensional analysis was performed using a me- 
thod outlined by Venezian (25) during the 1984 Texas A&M 
University dredging short course. 34 was chosen as the 
mee Of Length, $4°/v as the unit of time, and 0/34 as 
the unit of mass. The remaining seven variables expressed 


in these units give the following dimensional numbers: 
s/h, $,/H 3.1, Salde,, pl Gk = pr. elas, wine 
cegemeed ~1/3°? “d'-* %q’ ~50” s , dq & d 


By rearranging and combining terms: 


1 - Wy 
(S_,/D) (/S gu) v RP 


[ (S/h) (1/ (v/Sgu)*) (v?/8 478) 1* = u/-Yeh = Ny 
[(S4/D) (1/ (v/Sgu)*) (v2/8,43g)1* = u/VEd = Nop 


() (v2/8. 3g) (L/ (v/Sgu)*)1% = u/ V2848 = Neg 
[(1/ (v/S gu) 2) (v2/8 438) (Sq/Dgq) (/ (o, - 0) 17 
= U/ gDen (Co. - p)/e) = No 


Therefore, the functional relationship for scour Is: 


-. 


S 
—— = £(S4/Hy;3> Sg/D, Sq/D5q> Np» Np» Npp> Ngg- Nog): 








PD i 


PART IV 
EXPERIMENTAL EQUIPMENT 


The scour experiments were conducted in a 120-foot 
long, 3-foot deep and 2-foot wide two dimensional, glass- 
walled wave flume. The test area was located approximate- 
ly 40 feet from the wave generator. The entire test area 
measured 40 feet in length and was divided into several 
sections. An 8-foot long 5-inch deep test sediment bed 
was bordered on each end by an 8-foot long 10-inch high 
false bottom. At the leading edge of the test area, a 
16-foot long ramp was used to gradually bring the wave 
up to the new depth at the top of the false bottom. A 
schematic of the test area is shown in Figure (22). A 
1.5-inch aluminum pile was positioned along the center- 
line of the sediment bed. 

Currents were produced using a centrifugal pumping 
system. The flow rate was monitored using a Fischer and 
Porter electromagnetic flow meter. The pump discharge 
entered the wave flume just forward of the wave generator, 
and was diffused using a louvered grate. A Thermo-Systems 
model 1050 hot film anemometer was used to measure the 
vertical velocity distribution of the combined waves 
and currents. 

The random wave generator used was a Seeeun modular 


wave making system. This system consisted of a servo 
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control system amplifier, a Seasim programmable spectrum 
random signal generator, a Seasim autocompensating wave 
height gauge, and a paddle type sea wave simulation roll- 
ing seal modular wave maker. A capacitance-type wave gauge 
was used with a Hewlett-Packard model 17403A carrier pre- 
amplifier to measure the wave profiles. To determine the 
significant wave height from the wave profile, a wave 
digitizer was utilized. 

A wave absorber was positioned in the wave flume at 
the opposite end of the wave generator to reduce any re- 
frecaionsewhich might occur. 

Scour depth readings were made using a point gauge 


after each test. 





PART V 


EXPERIMENTAL PROCEDURES 


A total of 20 experiments were performed during this 
investigation. Four runs utilized currents only, and the 
remaining 16 runs used a combination of waves and cur- 
rents. All of the test runs used a 1.5-inch diameter pile 
placed in a bed of glass microbeads. The glass microbeads 
had a mean diameter (Ds 6) of 0.1 mm and a specific gravity 
of 2.45. For both the current only tests and the combined 
wave-current tests two water depths and two different cur- 
rents at each depth were used. Also used were four dif- 
ferent random wave spectra for each depth/current combina- 
tion. The data on the four random wave spectra are con- 
tained in Appendix A, and Table 1 contains a description 
of all tests conducted. 

At the commencement of each test, the sediment bed 
was leveled to a 5-inch uniform thickness. Water was 
added slowly to the tank so as not to disturb the sedi- 
ment. After the required water depth was attained, the 
proper flow rate through the tank was established. The 
water level in the tank was maintained by regulating a 10 
inch discharge valve located at the rear of the tank. 

Once the proper flow rate was attained and the required 
water depth held constant, the sediment bed was once again 


leveled. The signal generator, which had been previously 





TABLE 1 - SUMMARY OF TESTS 


Test Current Wave Spectrum 
IL U, 
Z U, 
3 U. 
4 U, 
6 Un D 
i U, E 
8 Un JL 
9 U, Ip 
10 U, le 
ii Uy vel 
eZ, Uy a2 
3 U, D 
14 Uy P 
IL U, ane 
16 U, J2 
1, U, D 
18 U, ie 
LS, U2 Jal: 
20 Uy JZ 
U, = .405 fps 
Us = .380 fps 
U, e650 fps 
U, = .485 fps 
D = Darbyshire 
P = Pierson - Moskowitz 
Jl = Jonswap (T = 8 sec.) 
J2 = Jonswap (T = 10 sec.) 
hy = 1.167 ft 
hy =e ol, £t 


2 as ae 
PP re 


SF oy 
Nm Nm 





programmed, was started generating the desired wave spec- 
trum. The maximum scour depth was recorded at various in- 
tervals throughout the run, the time increment increasing 
as the test continued. Each test was conducted until an 
equilibrium in the measured maximum scour depth was at- 
tained. Following the completion of the test, the wave 
generator, and the flow through the tank were stopped and 
the water was slowly drained from the tank in order to 
preserve the scour pattern formed. Once the tank was com- 
pletely drained, photographs were taken of the scour pat- 
terns. In addition, sediment elevations, taken in the form 
of a rectangular grid, were measured using a point gauge 
in order to determine the equilibrium scour depth. This 
depth was used in the data calculations presented later 


in this paper. 





PART VI 


PRESENTATION AND DISCUSSION OF RESULTS 


The experimental data obtained from this study are 


presented below in Table 2: 


TABLE 2 - EXPERIMENTAL DATA 


Run. Nelocrty (Eps) Scour Depth (ft) yt = 

1 U, = .405 101 

2 U, = .380 081 

3 U, = .650 167 

mA U, = 485 14 

5 U = .442 126 ys 

6 U = .394 ae 1656 
7 U = .393 13 1266 

3 U = .401 WDE 1715 

9 U = .386 | 09 1930 
10 U = .437 toy eer 
ili U = .459 152 1791 
12 U= .419 mie 2100 
13 U = .498 143 1724 
14 u = ol 149 ABD 
15 U = .504 174 “16eul 
16 U = .508 155 .1550 
17 U = .665 18 Migae 
18 U = .659 174 0967 
19 U = .660 176 1605 
20 U = .664 168 1093 


The results of the laboratory experiments are present- 
ed in graphical form in this section. These graphs depict 


the relationship between the dimensionless parameters 





developed earlier for scour as well as several dimensional 
plots. On various dimensionless plots, data from experi- 
ments conducted by Wang (22) and Armbrust (2) were used 

in combination with data obtained from the random wave- 
current interaction experiment. On several plots, curves 
were constructed to indicate a trend only, and these plots 
were labeled accordingly. Regression analysis was per- 
formed on the remaining plots in order to construct the 
appropriate curves. Logarithmic regression analysis was 
used on the semi-logarithmic plots, and power regression 
analysis was used on the full logarithmic plots. The Ween 
term adjacent to each curve represents the correlation 
coefficient for each curve. 

Figures (23) to (26) show the functional relationship 
that exists between the relative scour depth, Sq/h, and 
Sq/D and the sediment number for three distinct pile dia- 
meters and three different sediment sizes. It shows that 
scour depth increases with increasing wave-current velo- 
cities. Scour depth is not dependent on the sediment num- 
ber alone, but is also dependent on pile diameter (Figures 
23 and 25) and sediment size (Figures 24 and 26). 

Since the basic mechanism for scour hole development 
is the vortex system resulting from flow separation, and 
the separation point is dependent on the Reynolds number, 
the strength of the vortex system is a function of the 


pile Reynolds number. Because of the relationship between 
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the vortex system and the pile Reynolds number, scour depth 
should, in turn, be functionally related to the pile Rey- 
nolds number. This functional relationship, along with 
Shen's curve for scour due to a steady current, is de- 
picted in Figure (27). 

Figures (28) through (31) indicate the relationship 
between relative scour depth, S4/D and Sg/h, and the pile 
Reynolds number. It shows that scour depth increases with 
increasing wave-current velocities or pile diameter as in 
the case of the sediment number: it can be seen from the 
graphs that scour depth is not dependent on the pile 
Reynolds number alone, but is also dependent on the pile 
diameter (Figures 28 and 30) and the sediment size (Figures 
29 and 31). This agrees with the experimental work of 
Shen et al. (22); however, the dependence of the pile 
Reynolds number on sediment size differs from the conclu- 
sions reached by Palmer (18) in which he stated that the 
scour depth was independent of grain sizes from 0.12 mm 
eeO0-63 mm. 

By comparing Figures (23) to (26) and Figures (28) to 
(31) with figures found in Lee (13) which plotted the same 
variables for multiple piles, it is noted that the single 
pile condition requires a higher velocity to produce the 
same scour depth as the multiple pile case, assuming all | 


Sinem variables (i.e., Deg» h) are equal. 
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It was found that if the pile Reynolds number and 
the sediment number were multiplied together and plotted 
against the relative scour depth, SQ/h, all the points 
collapsed onto a single curve regardless of variance in 
pile diameter or mean sediment size. This result is 
shown in Figures (32), (33) and (34). The equation for 
each curve is also indicated on each figure. Due to the 
tendency of the data for random waves and currents and 
monochromatic waves and currents to collapse onto a sin- 
gle curve, it appears as if the irregularity of the waves 
is insignificant in terms of scour depth prediction. By 
using a dimensionless plot, as depicted in Figures (32)- 
(34), it appears as though predictions of scour depth 
could be made for situations in which the pile diameter 
or mean sediment size varied-without utilizing scaling 
factors. Before any conclusions can be made, field data 
are required to evaluate the need, if any, for scale fac- 
EOrS. 

Figures (35) through (37) show the relationship found 
for random wave-current interaction between the relative 
scour depth S q/Fy 73 and the sediment number, No, pile 
Reynolds number, Npp> and (N.) (Nap), respectively. The 
curves indicate that there is a rapid initial increase of 
the relative scour from the point of incipient motion up 


to a maximum relative scour depth. A further increase of 


N N 


5’ Npp oF (N.) (Nap) will cause a rapid decrease in the 
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Fig. 32. —- Relative scour depth, S4/h, vs. sediment 
number urmes pile Reynolds number for a constant 


pile diameter and various sediment sizes 
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relative scour depth, this decrease being independent of 
No: Mh ene (N_) (pp). These results are similar to those 
found by Wells and Sorenson (27) and presented earlier in 
this report in Figures (16) and (17). The relative scour 
depth is smaller in Wells and Sorenson's case, but this 
can be explained by the fact that they used wave action 
only and, therefore, their scour velocities were smaller 
in magnitude. The pile Reynolds number for the rapid de- 
crease of relative scour depth is approximately 7.5 x 10° 
sig! both cases, this rapid decrease in S/H ;3 being indepen- 
dent of the pile Reynolds number and the sediment size. 

A dimensional plot was constructed in Figure (38) 
showing the relationship of ultimate scour depth, Sq, ver- 
sus the combined wave-current velocity, U. It was found 
that in all cases the equilibrium scour hole was deeper 
than that associated with a steady boundary layer current 
with the same mean velocity. These results agree with 
predictions made by Niedoroda and Dalton (15). In gener- 
al, it can be stated that~the scour depth caused by wave- 
current interaction is approximately 10% greater than the 
scour depth due to a steady current alone. 

In Figure (39), scour depth was compared to Hy /3- 
Although the data are limited, it appears as though for 
current AINes ties greater than 0.45 feet per second, there 


is a tendency for scour depth to increase with increasing 


wave height, for velocities less than 0.45 feet per second, 
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scour depth decreases with increasing wave height and for 
velocities approximately equal to 0.45 feet per second, 
there is no definite trend at all. This result is most 
likely due to the fact that waves and currents interact 
differently as the magnitude of the current changes. Al- 
though a discussion was presented earlier in this paper 

on the mechanics of scour due to waves and currents, there 
is still a great deal of uncertainty in this area, and 


~Ueemner Lesearcen 1S) indicated. 


Scour Pattern Observations 


The rate of scour hole development by wave-current 
interaction was faster than that due to a steady current 
alone. A typical example of scour development is shown 
in Figure (40). As seen in this figure, scour development 
Hs very rapid initially, the scour rate then decreases 
until the equilibrium scour depth is reached. 

In all cases, the scour hole due to wave-current in- 
teraction was similar in shape to the scour hole formed by 
a steady current. This shape resembles an inverted cone. 
A typical example of a scour hole formed by a steady cur- 
rent is shown in Figure (41). Figures (42) to (45) depict 
the scour pattern formed by the same current used in Run 
4 (Figure 41), combined with four different random wave 
spectra. Photographs of the scour patterns formed by 
wave-current interactions for the remaining 12 runs can 


be found in Appendix B. As stated earlier, the equilibrium 
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Fige 42. - Scour pattern due to random waves and 
current — Run 13, 
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Fig. 43. - scour pattern 
Run i&4. 
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Fige 44. = Scour pattern due to random waves 
and current = Run 15. 
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Fig. 45. —- Scour pattern due to random waves 
and current — Run 16. 
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scour hole associated with the wave-current interaction 
was deeper when compared to the equilibrium scour hole 
formed by a steady current alone. The size and shape of 
the equilibrium scour hole associated with wave-current 
interaction was at times greater than, approximately equal 
to, or less than the size and shape of the scour hole due 
to a steady current alone. These results differ from the 
prediction made by Niedoroda and Dalton (15). They stated 
that the scour hole size and shape for combined waves and 
currents should not-be as great as that developed by a 
steady current alone. 

Scour development is illustrated by Figures (46) to 
(48). Figure (46) shows a typical pattern of scour develop- 
ment in its initial stages, approximately the first 2 to 
5 minutes of the test. Figure (47) and (48) are the typi- 
cal scour patterns that occur due to wave-current interac- 
tion 30 minutes and 2 to 3 hours after test commencement, 
respectively. 

The general scour pattern was one of localized scour 
around the pile. A vortex formed around the pile due to 
disturbed flow and suspended the sediment. The sediment 
was carried farther downstream under the combined effects 
of waves and currents than it was with a steady current 
alone. Initially, the sediment suspension appeared to 
become greater as the scour hole developed, possibily due 
to a strengthening of the vortex system surrounding the 


pile. 








Fig. 46. = Initial stages of scour development 
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Fig. 47. - Scour pattern 50 minutes after 
test commencement 
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48. = Scour pattern 2 to 5 hours after 
test commencement 


a= 








Because of the disturbance of flow caused by the pile, 
trailing ripples formed immediately behind the pile ex- 
tending in the direction of wave propagation as the test 
progressed. In nearly all cases, fairly uniform ripples 
formed along the entire sediment bed. The wave length of 
these ripples ranged from approximately 1.5 to 2 inches. 
Most of the ripples that covered the sediment bed initial- 
ly formed at the interface between the false bottom upstream 
of the pile, and the sediment bed. The ripples probably 
formed at this location due to a disturbance in flow, and 
once formed, the ripples progressed in the direction of 
wave propagation. 

In all cases, the maximum scour depth occurred down- 
stream of the pile, approximately 1/4 inch from the pile, 
and at an angle of approximately 40° to the left or right 
of an imaginary line drawn through the pile and parallel 


with the wave flume walls. 





PART VII 


CONCLUSIONS 


Scour development was very rapid initially. The scour 
rate then decreased until an equilibrium scour depth 
was reached. 

Scour development due to current and wave action oc- 
curred faster than scour due to a steady current alone. 
The equilibrium scour depth caused by wave-current 
interaction is approximately 10% greater than the scour 
depth caused by a steady current alone. 

The shape of the scour hole due to waves and currents 
and the shape of the hole due to current alone are 
similar. The shape resembles an inverted cone. 

For a constant pile diameter and sediment size, scour 
depth increases with increasing wave-current velocity. 
For a constant wave-current velocity and sediment siZe, 
scour depth increases a increase in pile diameter. 
For a current velocity greater than 0.45 feet/sec, 
scour depth increased with increasing wave height. 

For a current velocity less than 0.45 feet/sec, scour 
depth decreased with increasing wave height. 

A relationship exists between scour depth and sediment 
number (No). However, scour depth is not dependent on 
N. alone, but is also dependent on pile diameter and 


sediment siZe. 





9. A relationship also exists between scour depth and 
the pile Reynolds number (Npp) - Again, scour depth 
is not dependent on Np alone, but is also dependent 
on pile diameter, as well as sediment size. 

10. By plotting the relative scour depth (S3/h) against 
the product of the pile Reynolds number and sediment 
number , (Nap) (NS), all data points collapse onto a 
single curve regardless of variance in pile diameter 
or sediment size. 

ll. A relationship exists between the relative scour depth, 
(Sq/ Hy 3)> and Npp Ny and (No) (Nap) - 

12. In order to attain the same scour depth with a single 
pile as that achieved in a multiple pile case, the 
velocity must be greater assuming all other variables 


aGe 1dentical: 


4 


Although the amount of data is limited, due to the 
fact that all data could be collapsed onto one curve when 
plotting Sq/h versus (Npp) (No), (Figures 32 and 34), it 
appears as though the type of wave train used for model 
studies, either monochromatic or random, is insignificant 
in the prediction of scour depth about a single cylindri- 
cal pile. However, when attempting to predict the geometry 
of the scour hole, the irregularity of the wave train be- 
comes very significant. Lee (13) and Armbrust (2) stated 


that the pattern of scour for wave-current interaction 





is no longer conical in shape, but the pattern is irregular 
due to the wave particle orbital motion. The data in this 
study clearly indicate a conical scour pattern for random 
wave-current interaction. Although these results are dif- 
ferent from those of Lee and Armbrust, they agree with 
field data, as well as with the predictions of Niedoroda 
and Dalton (15). The fact that the geometric shape of 

the scour hole differs with the type of wave train used, 
may become important in the prediction of scour depth 
around multiple pile cases. 

Figures (35) and (36) show the relationship found for 
random wave-current interaction between the relative scour 
depth S4/Hy 3 and the sediment number, Ny» and pile Rey- 
nolds number, Nap» respectively. These results are simi- 
lar to those found by Wells and Sorenson (27) and present- 
ed in Figures (16) and (17). Because of similar trends 
in all figures and due to the fact that the pile Reynolds 
number for the rapid decrease of relative scour depth is 
approximately 7.5 x 10° in both cases, (Figures 36 and l/), 
it again appears as though the irregularity of the wave 
train is insignificant in the prediction of scour depth 
around a single cylindrical pile. 

In addition, Armbrust (2) found that there was an in- 
crease of scour depth with increasing wave height. As 
stated earlier, for irregular waves when the current ve- 


locity was greater than .45 ft/sec, scour depth increased 





L100 


with increasing wave height, and the scour depth decreased 
with increasing wave height for a current velocity less 
thanmeaoete/seck 

The conclusions resulting from this study were found 
for a flow in which the steady current was greater than 
the oscillatory flow produced by the waves. If the mag- 
nitude of the two flows in relation to one another change, 
the flow patterns will also change and the geometry of 
the scour hole may very likely differ from results found 
in this study. 

As mentioned in previous investigations, to attempt 
to predict scour depth for a prototype using the results 
of this study could lead to erroneous results. In order 
to accurately predict the prototype behavior, complete 
similitude is required between model and prototype, and 
field data is required to evaluate scale effects, should 


any exist. 


Suggestions EO meuliimathemeneseacch 


1. This study could be repeated under the same conditions, 
thereby omitting scale factors, if any exist, but using 
monochromatic waves to determine if the irregularity 
of the wave train is significant in scour studies. 

2. Further research can be conducted to investigate the 
variation of scour depth with varying significant wave 


height, Hy 3) as the current velocity varies. 





10) 1b 


It would be desirable to investigate why there is a 
rapid decrease in relative scour, S q/Hy 73> helmed pile 
Reynolds number, Nap» of approximately 7.5 x 10°. 

The study could be repeated using various sediment 
sizes and pile diameters to determine if scale effects 


are negligible when using a plot of S3/h versus (Np p) 
(N_). 


Ss 
This study could be repeated using multiple piles in 


order to compare the results with previous studies in- 
vestigating scour about multiple piles due to monochro- 


Matic wave-current interaction. 
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Random Wave Spectra Data 
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APPENDIX B 


Scour Patterns for Random Wave-Current Interaction 
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